Advances in nonequilibrium transport with long-range interactions 
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The effects of long-range interactions in quantum transport are still largely unexplored, mainly 
due to the difficulty of devising efficient embedding schemes. In this work we present a substan- 
tial progress in the interacting resonant level model by reducing the problem to the solution of 
Kadanoff-Baym-like equations with a correlated embedding self-energy. The method allows us to 
deal with short- and long-range interactions and is applicable from the transient to the steady-state 
regime. Furthermore, memory effects are consistently incorporated and the results are not plagued 
by negative densities or non-conservation of the electric charge. We employ the method to calculate 
densities and currents with long-range interactions appropriate to low-dimensional leads, and show 
the occurrence of a jamming effect which drastically reduces the screening time and suppresses the 
zero-bias conductance. None of these effects are captured by short-range dot-lead interactions. 
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Electron correlations have profound implications on 
the transport properties of nanoscale devices[l]. Lo- 
cal interactions within small molecules or quantum dots 
contacted to leads give rise to peculiar phenomena like 
Kondo effect [2] and Coulomb blockade [3], and have been 
the subject of several studies. Much less attention has 
been devoted instead to the nonlocal interactions respon- 
sible for interfacial screeening and polarization-induced 
renormalizations of the molecular levels. Recently, a 
short-range (SR) dot-lead interaction has been shown to 
cause a reduction of the quasiparticle gap due to the im- 
age charge effect [4-6]. In the interacting resonant level 
model (IRLM) the SR interaction is also at the origin of 
a negative differential conductance with a (interaction- 
dependent) power-law [7-10] as well as of an overall en- 
hancement of the off- resonance conductance [11, 12]. 

The theoretical progresses in dealing with SR dot-lead 
interactions are, unfortunately, not directly exportable 
to study long-range (LR) interactions, more appropri- 
ate for low-dimensional leads. The difficulty stems from 
the impossibility of combining many-body methods with 
embedding techniques, hence reducing the problem to 
the evaluation of the Green's function of a finite and 
interacting open system[13, 14]. Recently Elste and 
coworkers [15] approached the problem using the rate 
equations (RE) method in the IRLM with Luttinger liq- 
uid leads. The RE, however, are not reliable in the tran- 
sient regime and underestimate the steady-state polar- 
izability of the dot, as we will clearly show below. The 
fundamental questions which remain at present totally 
unanswered are therefore: What is the impact of a LR 
dot-lead interaction in the I-V curve? How does the 
screening time change from SR to LR interactions? 

In this Letter we consider the IRLM as the prototype 
model to address the above issues. We study the real- 
time evolution of the currrent and dot-density after the 



sudden switch-on of a bias voltage for both SR and LR 
dot-lead interactions. Our results indicate that LR inter- 
actions produce a jamming effect in the leads which (i) 
shortens the screening time and (ii) drastically suppresses 
the zero-bias conductance. 

The proposed methodology to conclude (i) and (ii) is 
based on a truncation of the equations of motion for 
dressed correlators. The procedure leads to Kadanoff- 
Baym-like equations with a correlated embedding self- 
energy which incorporates all interaction and memory 
effects. Our approach overcomes the negative proba- 
bility problem[16] of the RE and is, at the same time, 
charge-conserving. The final equations are exact in the 
uncontacted case as well as in the noninteracting case 
and several analytic results are obtained in the steady- 
state, including a Meir-Wingreen-like formula for the cur- 
rent. We benchmarked this formula against recent re- 
sults with SR interaction obtained using field theoreti- 
cal methods [7], DMRG[7, 12] and other renormalization 
group approach[9-ll], and found the same qualitative 
behavior. 

We consider the IRLM described by the Hamiltonian 
(in standard notation) 

+ J dxUix)pix)nd + Y,K€iO)d + h.c.], (1) 



with a = ±1 for R and L 



electrons, Ud = d^d and 
The dot-lead interaction 
U{x) in Eq. (1) can be either SR or LR. The system 
is driven out of equilibrium by the bias perturbation 
Hb J dxpa{x) = VaN^ (with the num- 

ber of electrons with chirality a) . For a non-perturbative 
treatement of the interaction we bosonize the fermion 
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operators[17, 18] 
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with boson field ^^(.t) = X^^^q ^^(^cw?^ ^'^'^^ — li.c.) — 
"^^J^" and an anticommuting Klein factor. In the 



mode expansion of the boson field A„ = 



with C 



"9 ~ v^' 

the length of the system and a a short-distance cutoff. 

The bosonizcd form of the electron density takes the form 

h.c.) + 
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N. 



-g-, and hence the bosonized Hamiltonian reads 
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(3) 



where Uq = J dxe"'"'U(x) and wo used U(x) = U{-x). 
Next we perform a Lang-Firsov transformation to (for- 
mally) eliminate the dot-lead coupling. The unitary op- 
erator W = e^^^°<! 2^^ii^lq-^''i)'^d transforms the orig- 
inal Hamiltonian into H = U^HU with 
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E 



h.c. 



(4) 



(from now on the sum will always be over q > 0). In 
the transformed Hamiltonian it appears the renormalized 
fermion field 



(5) 



evaluated in x 
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W, 



LLq 



- 0, with the effective interactions 

1 + Uq/ {2Trv) and WuLq = WLRq = 



Uq /{2Trv), and the renormalized energy level id = Sd + 

frtJZ'l^'jl^' basis the ground state of the 

isolated leads (i.e. for = 0) is the vacuum |0) of the bo- 
son operators baq- We can exploit this property to build 
the proper initial conditions by time propagation. We 
will consider the system initially uncontacted (T^ = 0), 
then switch on the contacts at time t = and let the cur- 
rent and dot-density relax. After relaxation, say at time 
to, we will bias the leads and study the screening dynam- 
ics from the transient to the steady state. This procedure 
simulates with high accuracy the so-called partition-free 
scheme[19, 20], as demonstrated in Refs. [21-23]. 

We define the dot Green's function on the Keldysh 
contour as 



1 



G{z,z') = -{T{d{z)dHz')}), 



(6) 



where T is the contour ordering, operators arc in the 
Heisenberg picture with respect to .ff -|- Hb (the bias 
perturbation does not change after the transformation), 
and the average is taken over the uncontacted ground 
state \0) (^\n), \n) being the state of the dot with single 
(n = 1) or zero (n = 0) occupancy. The Green's function 
obeys the equation of motion (EOM) 

{id,-id)G{z,z') = 6{z,z') + J2Ta{z)GM{z,z'), (7) 

a 

where Gaxiz,z') = j{T {fax{z)<l\z')}) is the dot-lead 
Green's function[24]. To close the EOM we derive Gax 
with respect to its first argument and find 

{idz + iavdx - Va{z)) Gax{z, z') 

= 7E<'^{[^/3*/V + l^-^-'/-] ^^)d\z)\). (8) 

The computation of the correlator in the r.h.s. of Eq. (8) 
is a formidable task. In order to proceed we approximate 
it by T*((/to/„, + /„,/y {z))pG{z,z'), where (...)p 

signifies that operators arc in the Heisenberg picture with 
respect to the uncontacted Hamiltonian. This approxi- 
mation is at the basis of our truncation scheme and be- 
comes exact in the non-interacting case as well as in the 
uncontacted case. Our approximation remains very accu- 
rate also for small since it correctly reproduces recent 
results with SR dot- lead interaction (sec below). 

To solve the EOM for Gax we define gaxaxi(z,z') = 

i(r {/a.(z)/L,(z')})p which satisfies the EOM 

{z,z') 

= S{z,z'){(u,fl, + fl,Ux) {z))p. (9) 

We can now perform a standard embedding and write 
the dot Green's function as the solution of 
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ed)G{z,z')- [ dzY,^c.{z,z)G{z,z') 



5{z,z'), 



(10) 

where Sa(z,z') = \Ta\'^gaOao{z, z') is the correlated 
embedding self-energy and the integral runs over the 
Keldysh contour. Using the Langreth rules[25] Eq. (10) 
is converted into a coupled system of Kadanoff-Baym 
equations (KBE) which we solve numerically. The real- 
time Keldysh components of E can be evaluated exactly 
using the bosonization method[17, 18] and read 
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with phase ipa{t) = Jq dtVa{t) and interaction dependent 
exponent 



Q{t) 



J2 — e-"9(e'''9* 
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FIG. 1. Time-dependent density for SR intraction using KBE 
(dashed), RE (solid) and MRE (dotted) for dot-energy id = 
(top panel) and Ed = 3 (bottom panel) and initial dot-density 
n = 0. The remaining parameters are Vl = Vr = 1, 7 = 0.2, 
U = 0.5, to = (partitioned scheme). Units are 10~^'i;/a for 
Va, 7, Sd, lOa/v for t and 2ttv for U. 



From solution of Eq. (10) we can easily calculate the 
dot-density from {nd{t)) = —iG^{t,t). Similarly, the cur- 
rent la Sit the interface between the dot and lead a can 
be calculate from 



^n{fUz)d{z)) +h.c. 
dz^a{z,z)G{z,z) +h.c. (13) 



In the steady-state regime G{t,t') depends only on the 
time difference t — t' and the current I = 1^(1 ^ 00) = 
—Init —5- 00) is given by a Meir-Wingreen-like formula 



da;S]>(c.)I]<(a.)-I]<(c.)S>(a.) 



2tt 



(14) 



Remarkably, the current cannot be written in terms of 
the difference between the leads Fermi functions despite 
the left and right contacts are the same. 

Our analysis starts by comparing the present approx- 
imation to the RE method, recently employed in a sim- 
ilar context [15]. In Fig. 1 we plot the time-dependent 
dot-density using the KBE, the RE and their Markovian 
version (MRE) for a SR interaction Uq — U. Both the 
KBE and RE densities exhibit oscillations with frequen- 
cies associated to charge-neutral excitations. As antici- 
pated, however, the RE suffer from the negative-density 
problem[16] (bottom panel). The MRE density is instead 
always non-negative but the lack of memory washes out 
the oscillations and the transient becomes a featureless 
exponential (top panel). The KBE density is superior 
also at the steady state. Both the RE and MRE predict 
a zero-temperature steady-state density either or 1 and 
hence severely underestimate the dot polarizability. We 
also verified that the KBE approach is charge-conserving 
since fulfills with high numerical accuracy the continuity 





0.0 




0.1 


— u = 


0.2 




0.5 




FIG. 2. 
7 = 0.2, 
for U. 



I-V curve for SR interaction. The parameters are 
id = 0. Units are lQ~^v/a for I, V, 7, id, and 2ttv 



equation d{nd{t))/dt — lR{t) + Ihit) at every time (not 
shown). 

We next calculate the steady-state current for a SR 
interaction and in the symmetric case {id — 0, Tl = 
Tn = T, Vl = -Vr = V > 0) recently considered by 
several authors[7, 9, 10]. In this case the integral in Eq. 
(14) can be performed analytically and, by defining the 
tunneling rate 7 — \T\'^/v, we have 



IsniV) ^ ( ^ 



f3~l 



1 



^r(/3) 
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Va 



1-/3 



-[ — ] np) 



with exponent /3(C/) — 1-1- 



(15) 

We notice that 



the functional form of /sr(^) is similar to the one de- 
rived in Ref. [9] within functional RG, although in the 
present case /3 is evaluated in a nonperturbative way. 
The above expression (plotted in Fig. 2) is also in ex- 
cellent agreement with the exact results of Ref. [7]. 
In particular it reproduces the universal ohmic behav- 
ior /sr(V^) — V/tt at small bias[26] (with ctq ~ I/tt the 
quantum of conductance), and the non- universal power- 
law decay /sr(T^) ~ V^~^ at large bias (the RE fail again 
here). The Authors of Ref. [7] observed numerically that 
the /3-exponent does not vary monotonically with J7, and 
for a special value U of the interaction reaches the maxi- 
mum value /3 — 1/2, for which the IRLM is exactly solv- 
able. Our formula, which is valid for all J7, is in fair good 
agreement with this result, and yields U — ttv. Note 
also that the steady-state current is symmetric around U 
since /3(C7 - 5U) = I3{U + SU). 

We can now present the most important numerical re- 
sults of the paper, i.e., the time-dependent current with 
LR interaction Uq — —Wlnlaq)"^. In this case the func- 
tion Q{t) as well as the integral in Eq. (14) must be eval- 
uated numerically. In Fig. 3 we display the I-V curve 
for several VF's. The behavior is qualitatively different 
from the SR case. In particular the zero-bias conduc- 
tance is strongly suppressed with increasing W. Due to 
the LR nature of the interaction the addition/removal 
of an electron to/from the dot induces a charge deple- 
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FIG. 3. I-V curve for LR interaction. Same parameters and 
units as in Fig. 2. W is in units of 2nv. 



tion/accumulation which extends smoothly deep inside 
the leads (jamming effect). For a current to flow the 
bias must be larger than the polarization energy of this 
particle-hole collective state. This picture also explains a 
common feature of the SR and LR I-V curves, i.e., the ex- 
istence of an optimal value of the interaction strenght for 
which the current has a maximum at fixed bias. Increas- 
ing the interaction from zero the electron density dimin- 
ishes close to the dot, thus enhancing the effective tun- 
neling rate (Coulomb deblocking). However, increasing 
the interaction further the particle-hole binding energy 
becomes larger than the charge-transfer energy Vl — Vr 
to move an electron from one lead to the other, and the 
current start decreasing. 

LR interactions have an impact also in the screening 
time. In Fig. 4 we plot the time-dependent currents for 
LR and SR interaction with same interaction strength 
W — U[27]. The LR current relaxes faster both in the 
partitioned scheme (contacts and bias switched on si- 
multaneously at t — 0) and partition-free scheme. The 
same behavior is observed for different values of W (not 
shown). The jamming effect of LR interactions is at the 
origin of the faster screening time. Electrons deep inside 
the leads suddenly respond to a change in the dot pop- 
ulation induced by the applied bias. Finally we observe 
that the steady-state value of the current is the same in 
both schemes. This agrees with the results of Refs. [20- 
23] according to which the memory of the initial state is 
washed out in the long-time limit. 

In conclusion we presented a comprehensive character- 
ization of the transport properties of the IRLM with LR 
interaction. We proposed an embedding scheme based 
on a suitable truncation of the EOM for the dressed 
fermion fields and derived KBE which we solved numer- 
ically and benchmarked against available exact results. 
The method was compared with recently proposed RE 
approaches, and found to be superior from the transient 
(no negative densities) to the steady-state regime (no se- 
vere underestimation of the dot polarizability) . LR in- 
teractions leave clear fingerprints in the time-dependent 
current as well as in the I-V curve, and we believe that 




Partitioned scheme 





LR (right) 




--- SR (right) 




LR (left) 




SR(left) 



Partition-free scheme 



— LR (right) 
-- SR(righO 



8 ^ 10 



FIG. 4. Time-dependent current for SR (dashed) and LR 
(solid) interaction for an initial density n = . The thick 
curves represent I Lit) while the thin curves /^(t). The pa- 
rameters are Vl = ~Vr = 1, 7 = 0.2, U = 0.2, Sd = 0. 
In the top panel to = (partitioned scheme) while in the bot- 
tom panel to = 60 (simulated partition- free scheme). Units 
are W~^v/a for Va, 7, id, Wa/v for t and 2'jtv for U and W. 



these features should survive when more sophisticated 
junctions (interacting multi-level resonant models) are 
considered. 
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